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CAS22 - Fortran Program for Fast Design and Analysis of Snock-Free 
Airfoil Cascades Using Fictitious-Gas Concept 

by Djordje S. Dulikravich* and Helmut Sobieczky** 


SUMMARY 

A user-oriented computer program, CAS22, was oeveloped that is applicable to 
aerodynamic analysis and transonic shock-free redesign of existing two-oimensiona l 
cascades of airfoils. CAS22 can be used in three separate modes of operation: (1) 

as an analysis code for full-potential, transonic, shocked or shock-free cascade 
flows, (2) as a design code for shock-free cascades that uses Sobieczky' s 
fictitious-gas concept, ana (3) as a shock-free design code followed automatically 
by the analysis in order to confirm that the newly obtained cascade shape provides 
for an entirely shock-free transonic flow field. CAS22 generates its own 
four-level, boundary-conforming grid of an 0 type. 

The mathematical model of the flow is a full-potential equation. Its 
artificially time-dependent form is solved in a fully conservative form by using a 
finite area technique; rotated, type-dependent upstream differencing; and successive 
line over relaxation. 

Isentropic shocks are captured by using a first-order artificial viscosity in 
fully conservative form. The shock-free design is performed by implementing 
Sobieczky's fictitious-gas concept of elliptic continuation from subsonic into 
supersonic flow domains. Recomputation inside each supersonic zone is performed by 
the method of characteristics in the rheograph plane by using isentropic gas 
relations. The new shock-free contour is determined from the condition that stream 
function is equal to zero on the airfoil surface. Besides being capable of 
converting existing cascade shapes with multiple snocked supersonic regions into 
shock-free cascades, CAS22 can also unchoke previously cnoked cascades and make them 
shock free. 


*Vi siting Research Scientist at NASA Lewis Research Center presently employed by tne 
Universities Space Research Association, Columbia, Maryland. 

**Research Scientist, DFVLR - Institute fur Theoretische Stroemungsmechamk, 
Goettingen, Federal Republic of Germany. 



INTRODUCTION 


Transonic flow is characterized by the coexistence of subsonic and supersonic 
local flow domains. These supersonic regions generally terminate with shocks, wnicn 
are the main generators of vorticity and entropy. As a consequence tne aerodynamic 
d rag force (waive 'drag )' 'and ^the; aerodifnami c noi se T e'\^'el i nahe'ase'; Vharply, and tn i s 
results in a rapid decay Of the' aerodynamic e^fficieh'ey-'df the “eaS When the 

shocks are strong enough, the shock - boundary layer interaction leads to boundary 
layer separation ano potential ly dangbrd^ iihs Wady' ( Vow and^ mechanical viorations. 

Besides these general features, transonic flow through tne cascades of airfoils 
is often characterized by the flow-choking phenomenon. Cnoking occurs wnen two 
neighboring airfoils in a cascade are connected by a single sonic line. Choking 
places an upper limit on the mass flow through a given cascade. It inevitably 
occurs in accelerated turbine- cascade flows with supersonic exit conditions. • 

FOr shock-free or snock-less flow the fluid decelerates from supersonic speed to 
subsonic speed hot discontinuOUsly (shocked flow) but smoothly ever- a: finite 
distance (iSentropic recompression.) . There are only a limited number of design 
techniques for obtaining shock-free cascade f lows. -. 'The method .of ; complex 
characteristics developed by Garabedian. (refs. - L and 2) and ,appl ied" to cascade - 
design, by Garabedian!and Korn (ref . 3) is probably the most matnenicitical ly rigorous 
and complex technique presently in use by industry (ref. 4). This method : was proven 
to work well in the design of high-performance isolated airfoils -and cascades^ ,lt . 
gives the designer , a freedom to prescribe velocity distribution along the entire 
(yet lunknown) contour iOf van- eirfoi 1 in the cascade. The result is an entirely ; 
shocktrfree. f low field and a cascade shape .compatible' with this f low, fie Iq., The mai-n 
drawbacks of this technique are that very often the input surface velocity , 
distribution is unrealisbicv and, the resulting cascade .shape does not have any , 
physical meaning. -This .technique cannot, be. applied tp. three dimensions, and for the. 
design of practically interesting shock-freeplanar cascades this. technique requires 
a sign.if icant ,aniount.-of user experience. An original and instructive technique for 
designing shock-free cascades of . airfoils .by uclog electrical analag was developed ' 
by: .Sobieezky (ref . 5) , - Because, of its, complex hardware requirements this unique 
technique hqs not found, wide acceptance in industry. 

A relatively simple method based bn the hodpgrapfi transformation was developed 
by Schmidt (ref. 6). It iS: fast and easy^to understand, but it has a. significant 
drawback associated with its stream function formulation. The resulting cascade 
shapes appear to have only sharp leading and trailing edges. As a consequence, the 
designer is forced to somewhat arbitrarily round these edges at the expense of the 
expected aerodynamic performance of the cascade. 

Yet another technique for designing shock-free cascade flows is presented in 
this report. It is based on the fictitious-gas concept of Sobieezky (ref. 7) and 
the full-potential, steady, transonic cascade analysis code of Dulikravich (ref. 

8). The present work has two main objectives. One is to eliminate ail shocks from 
the flow field. The other is to unchoke already choked cascades and to make tnem 
shock free. 

Both tasks will be achieved by partially reshaping ("shaving off") existing 
airfoil shapes in a cascade without changing tneir stagger angle, cnoro length, or 
gap-chord ratio. 

These goals seemed (ref. 9 and private communication with P. Garabedian, Courant 
Institute of Mathematics, New York University) to be accomplishable by prescribing a 
smooth, shock-free surface pressure distribution along a portion of a given airfoil 
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contour and then dete,riiiin.1r)g a. partially new ai rfoi 1 , shape that -i s cpnsi stent with ■ 
the prescribeo surface, pressure distribution. /because, ,pt ,tne high l,y fion linear ..." 
character of the.’transonVc 'flp'w this- design ,t,echnique genera,! ly does: not create, ari 
entirely shock-free flow field. ' 

It has been suggested in a number of standaro textbooks on gas dynamics 
(ref. 10) and accidentally observeo (private communication with J. McFaaen, ) that a 
number of shocks can remain in the flow field although each airfoil surface is shock 
free. These shocks have the unique property of terminating oefore reaching an 
airfoil surface. The conclusion is that an entirely shock-free flow fiela cannot oe 
obtained solely by prescribing shock-free surface pressure wi thput enforcing 
Shock-free conditions throughout the flow field. 

All these techniques were developed with a specific goal to create an entirely 
shock-free cascade flow field for a single set of global flow and geometric 
parameters. For this purpose a shock-tree surface pressure prescription tecnnique 
alone cannot be accepted as appropriate. 

Fictitious-Gas Concept 

For the specific purpose of designing shock-free transonic flow fields Sobieczky 
developed (ref. 7), elaborated (ref. 11), and successfully appliea (ref. 11!) trie 
idea of a fictitious-gas elliptic continuation. This truly ingenious concept is 
easy to understana and to apply, yet it is matnematica My and pnysically sound and 
correct. The method consists of three basic steps. The first step is an iterative 
determination of the shape of a sonic line oounaing a supersonic snocK-tree flow 
domain. This is attained by utilizing nonphysical (fictitious) expressions for trie 
fluid density and the speed of souno whenever the rlow is locally supersonic. 

The second step of the design procedure recomputes the flow fiela irisiae that 
supersonic domain by using now physical (isentropic) relations for trie fluid density 
and tne speed of sound. The third step then searches through tnis supersonic domain 
for the points. where stream function has a zero -value and thus determines a 
Shock-free airfoil shape that differs from the input airfoil contour only along the 
section wetted by the supersonic shock-free flow. Computer program CASZ2 unifies 
all three steps of the design procedure. 

Step 1 - The continuity equation for a steady, compressible flow normalized oy 
the critical values 



( 1 ) 


can be expressed in its canonical full-potential form (ref. 13) as 


p 

p* 







Here p is the density, u and v are the velocity vector components along tne x 
ana y coordinate axes correspondingly, M* is the local critical Mach numoer, s 
is the local streamline aligned coordinate and n is tne coordinate, locally 
orthogonal on the streamline. , When the flow is locally subsonic (M* < 1), 
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equation (2) is a nonlinear but elliptic. This means that the mass flow 
(pq/p*q*) will continue to increase with increasing of M* until the critical 
flow conditions are met (fig. 1). This can be concluded from the expression for 
mass conservation 


^ = M* ^ 
q* p* p* 


( 3 ) 


For an isentropic flow we can write 



(4) 


and 

a^/a^ = (p/p*)^ “ ^ 

where y is a specific heat ratio of a calorically perfect gas. Equation (3) has a 
maximum value of 1 when M* = 1. Once the flow field becomes locally supersonic 
(I'i* > 1), equation (2) changes its type and becomes nonlinear and hyperbolic. 

With the further increase of M*, mass flow will start to decrease. Supersonic 
flows generally revert to subsonic conditions discontinuously through a shock wave 
(or an isentropic discontinuity in the case of a full-potential flow model). To 
prevent this from ever happening, each time when the flow is locally supersonic 
(M* > 1) a fictitious-gas relation should be introduced so that the governing 
equation (2) continues to retain its elliptic character (ref. 14). 

This must be done in such a way that the mass flow never decreases after it 
reaches its physically maximum value of 1. Hence, the parabolicity condition 
(refs. 15 and 16) for any fictitious gas is 


That is, isentropic and fictitious-gas relations should coincide only along the 
sonic line where M* = 1. Recently, professor David Caugney of Cornell University 
pointed out an additional general constraint which every fictitious gas relation 
should satisfy. Namely, the slope of any fictitious density curve (fig. 2) must be 
less then the slope of the parabolic low density curve for any value of M* 
greater than one. 


dM* 

In addition, an analytic expression for the fictitious mass flow should have the same 
slope as the isentropic mass flow formula at M* = 1; that is. 
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( 8 ) 




d 

* 


* 



= 0 


(9) 


When M* = 1, it becomes 



From the general continuity criterion (refs. 16 and 17) one easily obtains a 
general formula for the fictitious speed of sound normalized with the critical speed of 
sound as 



For all practical purposes there is an infinity of choices for the analytic 
expression of fictitious density and fictitious speed of sound that all fulfill the 
conditions expressed in equations (6), (7), and (10). The main decisive factor 
should be the economy and speed of numerical evaluation of both fictitious 
quantities. When choosing a simple expression for pf/p* we must keep in mind 
that it should produce an equally simple expression for af/a». 

In the CAS22 program a simple parabolic-law expression was used 




( 12 ) 


■'which gives (ref. 16) 



1 - y'l + 4P(M* - 1) 
2P 


(13) 


and (from eq. (11) ) 
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(14) 




1 - yl + 4P(M* - 


2P 




y/1 + 4P(M* - 1) 


A single constant input parameter P has the purpose of making the fictitious gas 
more or less incompressible (fig. 2). This translates into the pnysical fact of how 
much longer and flatter the resulting supersonic shock-free flow domain will be than 
the original shocked supersonic zone (fig. 3). 

As already stated, the fictitious-gas relations are used during the iterative 
solution of the full-potential equation only at the grid points that indicate 
locally supersonic flow. This means that once the iterative procedure converges, 
the subsonic part of the flow field and the sonic lines will be correctly calculated 
from equations (1), (4) and (5). Hence lift and drag coefficients will be correctly 
determined, and the subsequent recomputation of tne supersonic domains will not 
change them. Note that these coefficients will correspond to the new shock-free 
cascade shape to be determined by design steps 2 and 3. These coefficients will 
differ from the coefficients of lift and drag tnat were obtained for the original 
shocked cascade by using isentropic relations only. 

The main drawback of the fictitious-gas concept as used in this work is tnat the 
user cannot prescribe tne specific form of the shock-free surface pressure 
distribution. The only thing that he can prescribe is the value of the constant P 
in equations (13) and (14) that will guarantee an entirely shock-free flow field. 

The resulting airfoil surface shock-free pressure distribution cannot be 
quantitatively influenced in a controlled fashion. 

It can be concluded now that the fictitious-gas technique snould be used for 
redesigning the existing transonic cascades and making them shock free. Ihe 
fictitious-gas concept is not limited in its application only to two-dimensional 
planar problems as are hodograph techniques; it has been successfully applied to 
both arbitrary two-dimensional (ref. 15) and tnree-dimensional (ref. 18) 
configurations. 

Step 2 . - After the iterative process for numerically solving the full-potential 
equation has converged, the entire flow field is searched for the location of sonic 
lines. CAS22 is capable of accurately determining the shapes of an arbitrary number 
of supersonic shock-free regions. It is not capable in its present version of 
detecting the shape of a choking sonic line obtained from the fictitious-gas 
computations. CAS22 can calculate flows where the leading edge can be covered 
completely by the supersonic flow, but trailing edge must be subsonic. Once the 
X*, y* coordinates of a sonic line enveloping tne particular supersonic zone are 
found oy interpolation, values of the velocity vector components u*/a* and v*/a* 
are determined together with the values of the velocity potential function (j>*. 

Values of the stream function \|i* along the sonic line are then determined by 
numerically integrating 

di|) * = oy* cos e* - dx* sin o* (lb) 


where 
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e* = tan ^ 


(16) 


is the local flow angle on the sonic line. Thus we obtained all the necessary data 
along the sonic line. These data will serve as initial data for the 
two-dimensional, method-of-characteristics marching integration (recomputation) of 
the full-potential equation inside the supersonic zone, where now only isentropic 
relations are used (ref. 7). The method of characteristics was used because it 
represents the most accurate, fastest, and simplest technique for integrating 
hyperbolic partial differential equations. The fact that we are integrating a 
quasi-linear hyperbolic equation is circumvented by rewriting the full-potential 
equation in its rheograph form (ref. 7). 



K(v)4> „ 





(17) 


Note that in the rheograph plane (local flow angle e versus Prandtl-Meyer 
function v(M*) the hyperbolic, quasi-linear, full-potential equation becomes 
linear. It is valid inside the triangular domain (fig. 4) bordered by the sonic 
line data (v* =0; «a £ £ ^b) ^nd the end characteristics 

Z = constant and n = constant, which are mutually orthogonal. Numerical 
integration of equation (17) can be easily performed inside the entire triangle ABC 
in the hodograph plane if the characteristic grid is uniform (fig. 4). This is the 
case only if the total flow deflection angle along the sonic line 

a-e* = 


is divided in a number of equal intervals (ref. 7). Values of x*, y*, 0*, and 

are then interpolated at these points along the sonic line by using cubic 
spline functions. Characteristics integration is performed in a line-by-line 
marching fashion from sonic-line initial data toward the point C of the triangular 
mathematical domain of dependence of the full-potential equation in a rheograph 
plane (fig. 4). 

Step 3 . - The mathematical domain of dependence should always be larger (fig. 5) 
than the physical, supersonic, shock-free zone that we are recomputing. With a very 
poor choice of parameter P in equations (13) and (14) it could happen that the 
mathematical domain of dependence will be smaller than the physical supersonic 
zone. Such a case is signaled automatically as an error, and the further execution 
of the CAS22 code is terminated. Yet another consequence attributable to the same 
cause is the occurrence of limit lines (ref. 19) inside the ABC mathematical domain 
of dependence. If limit lines do not protrude inside the physical supersonic 
domain, they are of no concern and the results of the computations are physically 
acceptable (ref. 14). If the limit lines occur inside the supersonic flow domain, 
the entire calculation should be repeated with a different (usually smaller) value 
of the input parameter P. 

Coordinates of the partially new shape of the airfoil in a cascade are 
determined from the condition that the stream function must have the value 
t|)s = 0 everywhere on the surface of the airfoil. The third step of the design 
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procedure searches through the entire triangle ABC for the pairs of points where 
i|) changes its sign. At these points <65 ®s S'"® interpolated, and the 

new airfoil surface coordinates Xg and between the points A and B are 

found by integrating the following expressions: 


cos sin ©g 


Sin © 


cos e. 


*^^S (p /p* ) (p/p*) 


(19) 


( 20 ) 


COMPUTER PROGRAM 
General Description 

Computer program CAS22 was developed with the intention to create a single program 
(fig. 6 ) that can operate either as an analysis code or as a shock-free design code 
for airfoil cascade flows. The analysis portion of CAS22 represents a modified 
version of the CAS2D analysis code for transonic cascades (ref. 8 ). In the design 
section of CAS22, subroutines CHASON and INVERS represent modified versions of their 
counterparts in the shock-free design code for isolated airfoils developed by 
Sobieczky (ref. 14). 

Computer program CAS22 consists (fig. 6 ) of the MAIN routine and subroutines. 
Subroutine MAIN reads the input data from unit 5 and generates x,y computational grid 
coordinates by calling subroutines CONMAP, SPLIF, INTPL, REMAP and XYINF. The grid 
(fig. 7 ) is generated very efficiently (refs. 8 and 20 ) by utilizing conformal 
mapping, elliptic polar coordinates, and nonorthogonal coordinate stretching and 
shearing transformations. CAS22 can generate a maximum of four consecutively refined 
grids. Subroutine CONMAP performs point-by-point conformal mapping from the physical 
x,y plane (fig. 7) into the X,Y computational rectangular domain (fig. 8 ). 
Subroutines SPLiF and INTPL represent simplified versions of the cubic spline-fitting 
and interpolation routines. Subroutine REMAP performs remapping from the X,Y 
computational rectangular domain back into the x,y plane. Subroutine XYINF modifies 
the positions of the grid points at upstream and downstream infinity and specifies x,y 
coordinates of imaginary grid points and periodic grid points. 

Subroutine ANALYS performs an iterative line overrelaxation of the continuity 
equation in its conservative form (eq. (1)). ANALYS actually generates coefficients 
of the three-diagonal correction matrix [C-jj] 

:c,j] rt.) . (Rj) + (Vj, (21) 


obtained from an artificially time-dependent form 


.((1 - t . 0 


( 22 ) 


(refs. 8 , 13, and 20) of the full-potential equation (2). The vector of residuals 
(Rjl is evaluated by applying a finite area technique (refs. 20, and 21) to 
equation (1). ANALYS uses isentropic relations only (eqs. (4) and (5)) and 
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captures all possible existing shocks (isentropic discontinuities) by performing a 
type-dependent, rotated finite differencing (refs. 8 , 13, 20, and 21) on the 
following artificially time-dependent form of equation ( 22 ): 

»(l - «*)<ss - "(l - * 2'=n^nt * '^t) ‘ ° 


( 23 ) 


Here H designates upstream differencing and E designates central differencing to 
be performed on the ^,xx» ^>YY» ^nd ^,xy derivatives constituting ^,55 
or (refs. 8 and 20). The artificial viscosity vector {Vj} 
represents the difference between upstream and centrally differenced <1 55 
derivatives in the first part of equation (23). 

As already mentioned, ANALYS captures isentropic shocks that differ in their 
strength and position from the Rankine-Hugoniot shocks. For the purpose of gaining 
a user's confidence in the capabilities of the CAS22 code to solve a full potential 
equation, it is advisable to compare the results of CAS22 with the exact results 
obtained iteratively from mass conservation (ref. 16) 

M* M* = M* ^ (24) 

a b a ^b 

A comparison of exact values for the isentropic and Rankine-Hugonoit jumps in Mach 
numbers across a normal shock is shown in figure 9 . 

Boundary and periodicity boundary conditions are applied explicitly in 
subroutine ANALYS and in subroutine BOUND. When CAS22 is to operate as a shock-free 
design code, it will automatically switch to subroutine XDES (fig. 6 ) instead of 
ANALYS. Routine XDES uses isentropic-gas relation at all the grid points where the 
flow is subsonic. At the points where M* > 1, subroutine XDES uses the 
fictitious-gas relations (eqs. (13) and (14)) for fluid density and the speed of 
sound. In such a way the full-potential equation is prevented from ever becoming 
locally hyperbolic, and the shock-free supersonic zones are created. Subroutine 
XDES thus represents step 1 of the shock-free design process described earlier. 

Step 2 is embodied in subroutine BUBBLE, which searches for the existence of shock 
free supersonic zones and determines the values of x*, y*, i|j*, and e* along 

their respective sonic lines. BUBBLE can perform this task on an arbitrary number 
of supersonic zones on a one-by-one basis starting at the trailing edge and moving 
clockwise around the airfoil. BUBBLE can treat even supersonic geometric leading 
edges, but it cannot treat choked shock-free flows. Step 3 of the design procedure 
involves subroutines CHASON and INVERS. CHASON performs characteristics integration 
of X, y, 0, and \|i inside the triangle ABC (figs. 4 and 5) in the rheograph 
plane by using isentropic relations. INVERS iteratively inverts the Prandtl-Meyer 
function (ref. 7). 

After the integration is over, CHASON searches for the values of Xg and 
yg where = 0 inside ABC. These are the new coordinates of the shock-free 
airfoil in the interval between points A and B (fig. 5). CHASON also searches for 
the existence of the limit lines inside the supersonic zones. After step 3 of the 
fictitious-gas, shock-free design procedure is finished, the CAS22 program proceeds 
with the closing subroutine CPMACH, or it automatically switches to its analysis 
mode. This switching capability is regulated by a single input parameter (explained 
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later). CPMACH calculates and prints surface Mach number, critical local Mach 
number, surface density, and the coefficient of pressure and plots (with the line 
printer) a rough plot of the Cp-distribution on the surface of the airfoil in a 
cascade. If the iterative calculation process is to be repeated on one or more 
additionally refined grids, subroutine MESH is called instead of CPMACH. MESH 
doubles the number of grid cells in each computational direction X and Y. MESH 
then interpolates values of tne potential function obtained on tne previous coarser 
grid onto the new grid by calling SPLIF and INTPL. In such a way an improved 
initial guess is created for the potential field calculation on tne new refined grid. 

INPUT 

The input data for CAS22 are identical to the input for the CAS2D program 
(ref. 8). The only difference is the addition of one computer card between cards 6 
and 7 in CAS2D input. This inserted card contains the following three additional 
input parameters: 

TITR4 maximum number of iterations to be performed on the final (fourth) grid by 
routine ANALYS or by routine XDES. In tne CAS2D program this value 
was predetermined to be one half of TITR3. 

SUPCEL number of equal angle increments to be generated along each sonic line 
between points A and B (figs. 4 and 5) before the characteristics 
integration is performed by CHASON. Tne maximum allowed value for 
SUPCEL is 299, although it is quite satisfactory to work with 
SUPCEL = iOO. 

FCTGAS fictitious-gas constant P to be used in equations (13) and (14). It is 
the most important input parameter in the CAS22 program because it 
also determines the operating mode of the program according to the 
following rules: 

FCTGAS > 0.11: CAS22 works as a design code 

FCTGAS > 0. out < 0.1: CAS22 worKs as an analysis code 

FCTGAS < 0.: CAS22 works as a design code followed automatically by the 

analysis mode in order to confirm that the newly obtained airfoil 
shape in a cascade is truly shock free. Actually the relation between 
FCTGAS and the constant P in equations (13) and (14) is 

P = 1 FCTGAS 1 

The user of CAS22 should be aware that he can to a certain degree 
influence the shape of the '-f^suT ting Cp-distribution. Very low 
(absolute) values of FCTGAS will barely smear out the original 
shocks. . Very high (absolute) values of FCTGAS will produce a highly 
"peaky" Cp-distribution (ref. 15) and will increase the cnance of a 
limit-line occurrence inside the supersonic flow domain. Generally 
acceptable values (absolute) of FCTGAS are 
5. < I FCTGAS 1 < 500. 

This new card in the input of CAS22 has tne same format (2x, 3 (7x, E12.5) as tne 
origin->l first few cards in the input of CAS2D (ref. 8). Examples of several inputs 
for the CAS22 program are shown in figures 11, 17, 20, 27, 32, 33, 34, 39. 

OUTPUT 

The output from CAS22 is almost identical with the output from the CAS2D program 
(ref. 8). Again, the present version of CAS22 does not have any need for tapes 
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although additional command for writing the results on a unit 14 can be easily 
incorporated. The difference between the output from CAS2D and CAS22 is that CAS22 
prints surface pressure, density, Mach number, etc., distribution only (ref. 8) 
after the solution has converged or the the maximum number of iterations has been 
reached on the final grid. CAS22 does not print out the Mach number chart (ref. 8), 
but it finishes the output with a rough printer-plot of the airfoil surface pressure 
coefficient calculated from the equation 

P - P CO 

^ 


with p_co = 1, o_o, = 1 and q_„ = 1. 

ASSUMPTIONS, LIMITATIONS, AND APPLICABILITY OF 
COMPUTER PROGRAM CAS22 

Computer program CAS22 is directly applicable for the analysis of transonic, 
steady, inviscid, and irrotational flow through two-dimensional static cascades of 
given airfoils. CAS22 is also directly applicable for the design (determination of 
coordinates of the points defining an airfoil shape) of shock-free airfoil shapes 
for transonic, steady, inviscid and irrotational cascade flow. It should be 
remembered that in both of these applications the identical program CAS22 is used. 
Only the inputs to CAS22 in these two operation modes vary. The difference in the 
inputs is confined to a single input parameter. 

In its analysis mode CAS22 can predict flow fields ranging in local speeds from 
incompressible through transonic, including the possible discontinuities (ref. 8) in 
the solution of the full-potential equation. These isentropic discontinuities do 
not represent physical shock waves because they do not satisfy Rankine-Hugoniot 
shock jump conditions (fig. 9). This suggests implicitly a basic assumption that 
all possible aerodynamic shocks must be weak (Mach number just ahead of the 
discontinuity should be Ma < 1.3) so that the entropy (actually vorticity) generated 
by them is negligible. The following assumption is that the boundary layer does not 
separate. Users should be aware that CAS22 does not calculate the boundary layer 
displacement thickness, the shock-boundary layer interaction, or the trai 1 ing-edge 
viscous effects. A separate computer program should therefore be added to CAS22 in 
the future that takes into account these phenomena. Cascades of closely spaced, 
highly cambered, thick turbine blades cannot be handled by the present CAS22 (refs. 

8 and 20) grid-generation technique. 

On the other hand, CAS22 can be successfully used for the analysis and shockless 
desiqn of isolated airfoils by giving the gap-chord input parameter a value of, say, 
PITCH = 30. 


RESULTS 

To demonstrate the capabilities of the CAS22 computer program in its design and 
flow analysis application, a series of numerical test cases (fig. 10) were conducted. 

For the purpose of simplifying input data preparation, coordinates of input 
airfoils were generated by using an analytic expression 

y = Ax + Bx*^ + ]/x(l - x) (D + Ex + F x(l - x) (26) 
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developed by Sobieczky (ref. 16 and DFVLR internal report, "Design Examples for 
Super-Critical Cascades", by H. Sobieczky and 0. S. Dulikravicn, to be published). 
Here A, B, C, D, E, and F are appropriate constants affecting the airfoil thickness, 
the curvature of the leading and trailing edges, and the mean camber line. 

Numerical flow calculations for each test case were performed on four successively 
refined, boundary-conforming, geometrically periodic grids of an 0 type. 

After 200 iteration cycles were performed on the first grid, the potential field 
was interpolatea on the second grid, on which an additional 200 iteration sweeps 
were performed. The potential function was then interpolated on the third grid, 
where 100 iterations were performed. Finally tnis potential field was interpolated 
on the fourth grid, where it served as an initial guess for the concluding 100 
iterations (TITR4 = 100.), after which the potential field was considered converged 
for the engineering applications. The grid sequence used in both the analysis and 
design modes was the following: first grid, 24 x 6 grid cells, second grid, 48 x 12 

grid cells, third grid, 96 x 24 grid cells, and fourth grid, 192 x 48 grid cells. 

The number of refined computational grid cells along each sonic line (fig. b) was 
SUPCEL = 98. Wake and boundary layer displacement thickness were not taken into 
account. The specific heat ratio used in the isentropic relations was y = 1*4. 

All calculations were performed on an IBM 3033 computer. 

Test case 1 is an airfoil cascade with tne following specifications: stagger 

angle e, -50°; free-stream angle at^upstream infinity a\, 65°; free-stream 
angle at downstream infinity 42“; Mach number at upstream infinity Mj^, 

0.80; gap-chord ratio g/c, 0.8. 

The entire input for test case 1 is shown in figure 11. An analysis mode of the 
CAS22 program with FCTGAS = 0. shows (figs. 12 and 13) that this highly staggered 
cascade flow has a single supersonic zone terminated by a strong shock. A design 
mode of CAS22 (with FCTGAS = 25.) "shaves off" this original airfoil cascade and 
creates a new shock-free (figs. 14 and 15) cascade flow field. It should be pointed 
out again that the actual change in the airfoil snape is performed only along tne 
section wetted by the supersonic flow. From figure 16 it can be concluded that the 
changes are infinitesimally small and reflect the extreme sensitivity of the 
transonic cascade flow to any changes in global flow conditions or in the cascade 
airfoil shape. Figure 16 also shows that the entire flow field is shock free. 

For the industrial application of shock-free cascades it is very important to 
know how tney perform at off-design conditions. I'o provide an illustrative answer 
to this question, coordinates of the newly obtained shock-free airfoil were 
substituted in the input data (fig. 17) set with an off-oesign, free-stream Mach 
number at upstream infinity M]^ = 0.78. Witn this input (fig. 17) an analysis 
mode (FCTGAS =0.) of CAS22 shows that the shock-free cascade reverts to a shocked 
cascade at off-design conditions. Tne question is now; Is this shock stronger or 
weaker then the shock that exists in the original cascade at the same off-design 
conditions? The answer is presented in figures 18 and 19. It is apparent that the 
shock-free airfoil cascade performs better then the original shocked cascade even at 
off-design conditions. The shock-free cascades will have shocks at off-design 
conditions, but these shocks will be weaker and the losses due to shock - boundary 
layer separation ano shock wave drop will be less. In adoition, tneir supersonic 
zones will be flatter, thus delaying the occurrence of choking flow conditions. 

Test case 2 (fig. 10) has the purpose of showing that CAS22 is capable of 
converting relatively thick airfoil shapes into shock-free snapes. Input data for 
test case 2 are shown in figure 20. With this input CAS22 reveals (figs. 21 and 22) 
that the original cascade has two supersonic zones terminated by strong shocks. 
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This cascade was then converted into a shock-free cascade (figs. 23 and 24) tnat nas 
a remarkably thin and elongated supersonic zone on the suction surface (fig. 2b). 

To show the improvement in aerodynamic efficiency of the cascade over a range of 
free-stream inlet Mach numbers, this new shocK-free cascade was analysed (fig. 26) 
at a number of off-design Mach numbers. One of the input data sets for these 
analyses is shown in figure 27. For weak shocks the shock wave drag losses are 
approximately proportional to the third power of the drop in Mach number aM across 
the shock. It is clear from figure 26 that the shock free cascade designed oy using 
the fictitious-gas concept offers better performance over a range of inlet Mach 
numbers although the snock-free design was performed only for a single value of 
Ml = 0.72. 

Test case 3 has the purpose of showing that CAS22 is capable of both analysis 
(figs. 28 and 2y) and a fictitious gas, shock-free redesign (figs. 30 and 31) of 
high-solidity (g/c = 0.6666), moderate-stagger (e = -25°) airfoil cascades. Note 
that the original cascade has a supersonic zone starting very close to the leading 
edge. Input data for the redesign of the original cascade in test case 3 are snown 
in figure 32. Input data for the analysis of tne newly obtained snock-free cascade 
in test case 3 are shown in figure 33. The value of FCTGAS used in this test case 
(fig. 10) was somewhat lower, because the higher values would cause very peaky 
surface pressure distribution and possibly limit lines inside the supersonic part of 
the flow field. 

Test case 4 has the purpose of showing in a simple way the capability of CAS22 
to unchoke already choked, shocked cascade flows and to make the unchoked flow 
fields entirely shock free. A nonstaggered cascade of NASA 0018 airfoils was 
analysed with Mi = 0.582. Tne flow is obviously choked (figs. 34 and 35) because 
even one-dimensional analysis predicts (ref. 22) that this flow will choke if 
Ml > 0.577. A design mode of CAS22 with P = 500. (eqs. (13) and (14)) converted 
this flow field into the shock-free unchoked cascade flow shown in figures 36 to 

38. Coordinates of the new shock-free NACA 0018 airfoil cascade are shown in figure 

39. 

The computer time required for each test run varied between 290 CPU seconds and 
330 CPU seconds on an IBM 3033 computer. Total core memory required was on the 
average 500K. 


CONCLUSIONS, SUGGESTIONS, AND 
PERSPECTIVES 

The fictitious-gas concept of Sobieczky was successfully applied in a design 
procedure for shock-free, two-dimensional airfoil cascades. Computer programs like 
CAS22 can be relatively simply developed from the existing full-potential analysis 
codes for compressible (not necessarily supersonic) flow fields. CAS22 represents a 
fast computer program ready for use by aerodynamic! sts with practically no knowledge 
of numerical methods and computer programming. A very concise and simple to 
understand input format is designed specifically for such users. The entire code is 
well documented by the comment cards, and possible changes or modifications can be 
easily added or inserted. In its present version CAS22 does not have a need for any 
tapes or disks. The entire output, containing coordinates of the new shock-free 
airfoil surface and the corresponding values of coefficient of pressure, density, 
Mach number, and critical Mach number, is printed on paper. If one desires to save 
them permanently for the purpose of a computer graphics display, WRITE(14) cards in 
the routines MAIN and BUBBLE should be reactivated. Each of them has an adjacent 
self-explanatory comment card. It should be clearly pointed out tnat the locations 
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of airfoil leading-edge and trai 1 ing-edge stagnation points are not fixed in CAS22 
but will change during the course of iteration. They will finally be forced to 
settle at the locations on the airfoil surface that are compatible with the input 
data parameters, ai, ap, Mj, 6, and g/c. For closely spaced olaaes at 
transonic speeds, the flow field at the trailing edge is mostly influenced by the 
value of ap. Nonrealistic values of any of these input parameters and especially 
of ap will force the trai 1 ing-edge stagnation point to move away from the 
geometric trai ling-edge location. This will result in excessively high supersonic 
Hach numbers at the trailing eage as the flow is forced to go around it. As a 
result the CAS22 iterative procedure will quickly diverge and the program execution 
will terminate. The iterative calculation performeo by CAS22 sometimes terminates 
if the user attempts to calculate choked flows, although it does not always happen. 
Termination or simply oscillatory iteration can also occur if one tries to redesign 
a shocked cascade with the Mach number just ahead of the shock higher then 
approximately = 1.35. 

If the user is not satisfied with the resulting shock-free surface pressure 
distribution, he can try two things; either repeat the run with a different value 
of the input parameter FCTGAS, or repeat the run with a slightly changed input 
airfoil shape. The second option gives a user an almost endless possibility to 
generate significantly different shock-free surface pressure distributions by simply 
adding shallow "bumps" to the original airfoil surface so tnat CAS22 will have more 
freedom to perform its "shaving" procedure. 

CAS22 does not take into account boundary layer displacement thickness although 
It can include the global wake blockage effect (ref. 8) via the input parameter AR. 

Any shock free cascade flow designed by CAS22 will be theoretically shock free only 
for a single set of input parameters. If any one of these parameters is changed, 
shocks will occur. Use of CAS22 in connection with an optimization program should 
generate airfoil shapes that offer the best aerodynamic efficiency (weakest shocks) 
over a specified range of, say, inlet Mach numbers. 

Use of existing analysis codes to clarify that the shock-free flow field is 
actually shock free as predicted by CAS22 can give slightly varying results 
depending on the accuracy of the numerical scheme used (ref. 23) and the grid 
resolution used in the particular analysis code. Figure 40 shows a comparison 
(ref. 16) between the results obtained from a design mode of CAS22 and the results 
of analysis obtained by the analysis mode of CAS22 and by an Euler equation solver 
(provided by Dr. David Ives of Pratt and Whitney Aircraft Company aerodynamics research 
staff) developed by Ni (ref. 24). CAS22 analysis version gives slightly sharper 
recompression due to the first order artificial viscosity used in ANALYS. 

Extension of the applications of trie fictitious-gas technique to 
quasi-three-dimensional stator shock-free blade design was recently accomplished 
(ref. 25). The final goal is a fully three-dimensional, snock-free rotor blade and 
nub design capability. 
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Figure 3. - Example of the sonic line shapes compatible with the 
isentropic shocked flow, parabolic shock-free limit, and a shock-free 
flow obtained by using fictitious-gas relations 








Figure 5o - Characteristic grid in physical and 
rheograph plane with shape of new (shock free) 
surface 
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Figure 6. - Global flow chart of CAS22 program 
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Figure To - Computational grid in physical 
(x,y) plane 
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Figure 8o - Periodic flow field in computational (X,Y) plane 
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Figure 9, - Rankine-Hugoniot 
(Mb ) and Isentropic (Mb , M^ 
shock jump relations 


25 




ANALYSIS AND DESIGN INPUT PARAMETERS 

OFF-DESIGN 

TEST 

CASE 

Ml 

(FMACH) 

3° 

(TWIST) 

°1° 

(ALPHAl) 

0^0 

(ALPHA2) 

G/C 

(PITCH) 


Ml 

(FMACH) 

■■ 

0.80 ■ 

-50. 

65. 

42. 

0.8 

25. 

0.78 

2 

.72 

-27. 

37. 

16. 

1. 

30. 

.70/. 75 

3 

.75 

-25. 

5 

13.5 

.6666 

2 0. 


■ 

.582 

0 

1 

0 

0 

1. 

500. 



Figure 10 o - Test cases for CAS22 program 
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11. - Test case 1 input data set (design mode) for original 
cascade 
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Figure 15, - Test case 1 Mach number 
distribution for designed shock-free 
cascade 
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Figure 17. - Test case 1 input data set (analysis mode) for shock- 
free airfoil cascade 
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20. - Test case 2 input data set (design mode) for original 


cascade 
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Figure 21. - Test case 2 coefficient- 
of-pressure distribution on original 
cascade surface 
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Figure 27. - 
free airfoil 


Test case 2 input data set (analysis mode) for shock- 
cascade 


43 




Figure 28, - Test case 3 coefficient 
of pressure distribution on original 
cascade surface 
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Figure 29. - Test case 3 Mach number distribution on 
original cascade surface 


45 












TEST CASE 

NO. 3 : ORIGINAL CASCADE 
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Figure 32. - Test case 3 input data set (design mode) for 
original airfoil cascade 
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TEST CASE 
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Figure 33. - Test case 3 input data set for shock-free cascade 
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TEST CASE 


XCELL - 
ALPHA1= 
PITCH = 
FMACH = 
TITRl = 
TITR^ = 
POINTS= 


0 00 0 00 
9911^9 
9SS699 
986111 
9832^6 
979970 
976153 
971670 
966^02 
960239 
953082 
9A^539 
935«»30 
92A789 
912861 
0.899603 
0 . 88^988 
0.869002 
0 .8516^^t 
0.832929 
0.812885 
0.791559 
0 . 768990 
795260 
720995 
699635 
667 930 
690991 
612286 
5835C9 
559982 
525097 
995579 
966051 
936667 
907559 
378852 
3507 07 
323218 
29651 3 
27 07 03 
295888 
222159 
1 99595 
178263 
158219 


NO. 9 : NACA0018 
29000D+02YCELL 
OOOOOD+02TWIST 
OOOOOD-^OOROl 
58200D+00RLX 
OOOOOD+02TITR2 
00000D+02SUPCEL 
69000D+02GAnMA 
0.000000 
0.001873 
0.002388 
0.002930 
0.003528 
0.009210 
0.005000 
0.005929 
0.007002 
0.008256 
0.009699 
0.011397 
0.013207 
0.015286 
0.017586 
0.020109 
0.022839 
0.025768 
0.028891 
0.032189 
0.035691 
0 .039225 
0 . 092917 
0 . 096688 
0.050510 
0.059351 
0.058176 
0.061951 
0.065638 
0.069197 
0.072589 
0.075772 
0 .078707 
0 . 081353 
0.083672 
0 . 085626 
0 . 0S71S5 
0 . 088319 
0.089005 
0 . 089229 
0 . 088980 
0 . 088256 
0.087065 
0.085918 
0 . 083339 
0 . 080855 


ORIGINAL CASCADE 
0.05000D+02PMESH 
O.OOOOOD+02ALPHA2 
0.01580D+OOR02 
1.70OOOD+OOAR 
2.00000D+02TITR3 
1.98000D+02FCTGAS 
1.90000D+00CONVER 


9.00000D+00 

0 . 000000-^02 

0.00050D+00 

l.OOOOOD+00 

l.OOOOOD+02 

5.00000D+02 

l.OOOOOD-08 


0 . 139509 
0.122197 
0.106161 
0 .091596 
0 .075288 
0 . 066359 
0 .055718 
0.096310 
0.038067 
0.030909 
0 .029797 
0.019979 
0.019995 
0.011178 
0 .007902 
0.005038 
0 .002999 
0.000000 


-0 . 077999 
-0.079813 
-0.071390 
-0 . 067626 
-0 . 063720 
-0 . 059672 
-0.055529 
-0.051337 
-0 . 097138 
-0.092967 
-0.038851 
-0,039800 
-0 .030805 
-0 .026819 
-0 .022733 
-0.018302 
-0.012886 
0.000000 


Figure 34. - Test case 4 input data set 
(design mode) for choked cascade of 
NACA 0018 airfoils 
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Figure 38. - Test case 4 sonic line shapes foi 
shock-free unchoked cascade of originally 
0018 airfoils 
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TEST CASE NO.^ : NACA0018 SHDCKLES5 CASCADE 
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Figure 39. - Test case 4 input data set (analysis mode) for a 
shock-free unchoked cascade of an originally choked cascade 
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Coefficient of pressure. 


Shockless design mode of CAS22 code 
Analysis mode of CAS22 code 
Analysis with an Euler solver 



Figure 40.- Comparison of results of exact shock-free design with results of analysis 
mode of CAS22 code and Euler equation solver (ref. 27,28) 
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